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A generalization to the theory of massive gravity is presented which includes three dynamical 
metrics. It is shown that at the linear level, the theory predicts a massless spin-2 field which is 
decoupled from the other two gravitons which are massive and interacting. In this regime the 
matter should naturally couple to massless gravitons which introduce a preferred metric that is 
the average of the primary metrics. The cosmological solution of the theory shows the de-Sitter 
behavior with a function of mass as its cosmological constant. Surprisingly, it lacks any non-trivial 
solution when one of the metrics is taken to be Minkowskian and seems to enhance the predictions 
which suggest that there is no homogeneous, isotropic and flat solution for the standard massive 
cosmology. 
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I. INTRODUCTION AND SUMMARY 

The general theory of relativity (GR) has been reigning supreme since its introduction by Einstein nearly a hundred 
years ago. The underlying principles upon which GR is founded, namely the principles of equivalence and general 
covariance, have given the theory, in spite of its inevitable shortcomings, an unprecedented breath and depth, enabling 
it to predict, solve and answer questions which seemed beyond the possibilities of any scientific endeavor over the 
past decades. That said, new observations and discoveries in the recent past seem to be beyond the predictive powers 
of the standard GR. The discovery of the accelerated expansion of the universe relating to dark energy and galaxy 
rotations curves relating to dark matter are some of the examples that necessitate modifications to the standard GR. 
Among many proposals to modify GR, Massive Gravity (MG), whose roots go back to the work of Fierz and Pauli 
[1] in 1939 has been gaining momentum in recent years after being brought back from oblivion in the 1970's where 
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renewed interest in quantum field theory became pervasive. It is an effective field theory of gravity which, as the 
name suggests, is an extension of gravity with a massive spin-2 graviton. Its linear version coupled to a soiirce was 
studied by van Dam, Veltman and Zakharov [2] who discovered that it makes predictions different from the linear 
GR even in the limit of the zero mass. This is known today as the vDVZ discontinuity. This problem is traced back 
to the degrees of freedom entering the theory, 5 for a massive and 2 for a massless spin-2 graviton. Attempts to 
remedy this shortcoming is suggested in [3]. However, Boulware and Deser [4] demonstrated that the theory suffers 
from the existence of ghosts, when non-linear self-interaction terms added. A fix was suggested in [5] where it was 
demonstrated that the resulting theory is unitary and ghost free. A comprehensive and well written review on this 
subject can be found in [6] . 



Recently, it has been proposed in [7] that massive gravity in a ghost free form can be written in a bimetric 
language. Both of these topics separately have had their supporters and followers. Bimetric theories are theories with 
two dynamical spin-2 fields. However, it has been shown that there arc no consistent theories of interacting massless 
spin-2 fields. In general these bimetric theories suffer form having Boulware-Deser ghosts, and so arc unstable [4]. 
The bimetric theory was first proposed by N. Rosen in [8] to give a tensor character to quantities like gravitational 
energy- momentum density pseudo-tensor. In this theory the tensor g^i, is regarded as a gravitational field that has 
no direct connection to geometry and so geometrization of gravity is given up [9]. Subsequently a large amount of 
work has been done to address cosmological problems in this context [10, 11]. 

Contrary to N. Rosen in [8] and [9], some of these works have kept the geometrical viewpoint of General Relativity. 
Milgrom has used two metrics for constructing relativistic formulation of MOND gravity [11] which naturally leads 
to MOND and dark energy effects. In this BIMOND theory, matter lives in the space-time described by one of the 
metrics which couples to another metric through an interaction term. Drummond in [10] introduced two vierbein 
bundles into the space-time manifold so that each bundle supports its own metric. One of this metrics is associated 
with matter and the other with gravity. This theory is a kind of modified theory of gravity that has the flexibility 
to permit the introduction of a length scale to explain the observed rotation curves of galaxies. On the other hand, 
in the last two years, the theory of massive gravity has opened a new era of research after the work of de Rham, 
Gabadadze and ToUey [12] where a new interaction term for a non-linear massive gravity was constructed, leading to 
a theory which becomes ghost free. The cosmology of such a model is considered and shown to be devoid of the flat 
FRW solutions [13]. However it has recently been shown that the open FRW solution does exist [14]. 

In ordinary formalism of massive theories of gravity an additional reference metric f^i, is required. It has been 
shown recently by S. F. Hassan and R. A. Rosen [7] that if we use the interaction term constructed in massive gravity 
theories, resulting in a ghost free massive theory, one will arrive at a bimetric theory of gravity which is free of the 
Boulware-Deser ghost. This model is important in the context of bimetric theories due to the form of deflnition of 
the potential term. The cosmology of such a bimetric theory is studied in [15, 16]. 

In this work we generalize the massive gravity scenarios to include three dynamical metrics. It will be clear that this 
trimetric formalism is very similar to A^-metric models, though it has new predictive powers as in bimetric models. 
Naturally, by introducing N metrics. A'' gravitons are expected in its linear regime. This fact is similar to the Kaluza- 
Klein modes. In this context the existence of A^ gravitons and the interaction between them are considered to address 
an effective field theory for gravity [17, 18]. In [17] there is an interesting interpretation of massive gravitons using 
the Fierz-Pauli mass in that the authors assume to have A^ sites which arc related by interaction through the mass 
terms. This interpretation can also be used in our assumption of N metrics, not in the Fierz-Pauli format but in the 
new ghost-free massive gravity formalism. 

As mentioned above, at the linear level, we arrive at a theory which describes one massless spin-2 field which is 
decoupled from two interacting massive spin-2 fields. In addition, the cosmological solution of this model will be 
presented. This model, similar to the very recent results in bimetric cosmology [15, 16], has some de-Sitter solution 
with a cosmological constant which is a function of the masses. It is interesting that in this model which is more 
extensive than ordinary massive cosmology [13], there is no non-trivial cosmological solution if one assumes one of the 
metrics to be Minkowskian. So it seems that the result of [13] can be generalized to: "there is no non-trivial solution 
for massive cosmology (even in its multi-metric representation) by assuming one Minkowski metric" . 



We start with a review of the non- linear massive gravity action presented in [12] . The non-linear ghost-free action 
which reduces to the Fierz-Pauli action at the linear limit is 



II. THE MODEL 
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where the square root matrix is defined as yfAyTK = A for a general metric A. The efc(-\/ g~^f) are 5 elementary 
polynomials of the eigenvalues A„ of the matrix y/g^^f. They are explicitly written as 



eoiVF^) = 1 



62 



1<J 



i<j<k 
4 

eiiVg^) = 11^ = detVg^ (2) 

i=l 

This implies that the highest order term in the interaction part of the action is simply a determinant of the background 
metric ffi^. We also note that one can write the polynomials &k{^/ g~^ f) in terms of trace of the matrix \J g~^ j ■ One 
observes that the potential term in the action is symmetric under the transformation [7] 

/ o 5, ^„ o /34-„, (3) 

so one may consider the above interaction term to be the interaction for the metric g^i,^ with the background metric 
j^y or the interaction for the metric /^,^. These observations led S. F. Hassan and R. A. Rosen [7] to add a kinetic 
term for the metric /^j^ to the action and converted it to an action for a bimetric massive theory which is found to 
be ghost free. In the present work we generalize the work in [7] to become a non-linear ghost free action describing a 
trimetric model. The action then reads 

S = -M^ j d^xy/^Rg -M] j d'^xyf^Rf -Ml j d-^xV^Rh 
f 4 

71=0 

r 4 

+ 2rniMg\ / d^Xy/^ ^ l^e^W g-^h) 

m=0 

f 4 

+ 2mlM% / fi4a;x/^Ea.e,(/f^) (4) 

•' s=0 

where Eg, Rj and Rh are the Ricci scalars constructed by three metrics g^i^, and /i^,^ respectively. We also have 
introduced three different Planck masses for each metric. We have for each interaction term, an effective Planck mass 
constructed by metrics included in the interaction 

1 1 1 
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III. LINEAR THEORY 



As is well known, bimetric theories in general, describe one massless and one massive spin-2 fields. In this section 
we will see that the above trimetric theory describes one massless spin-2 field and two massive spin-2 fields which in 
general interact with each other. We assume for simplicity, that the interaction terms are all given by the minimal 
model introduced in [7, 19] 



/3o = 3, A = -l, /32 = 0, /33 = 0, 134 = 1, 



(6) 
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and also for 7„ and a„. In this model the interaction terms can be written in terms of trace and determinant 



d^xy'—g 



(^3-^r^/^ + c^e^A^), 



and similarly for the other interactions. 

Now we expand the metrics around the same fixed background[20] 

To second order in perturbations, this trimetric minimal model reduces to the action 

- \miMif j d'x - r.y - {r, - nY 

- \mlMl^j d^x [(fft - KY {r, - 

- \mlM% j d'x [(/> - h^^Y - {f''^ - '^'^^Y 



(7) 



(8) 



(9) 



where is the Einstein-Hilbert kinetic term and the background metric g^i, is responsible to raise and lower the 

indices. Now by using the following transformations 



- ^ {mX, + M}j^, + M2 V) (10) 

^HU = Qiiu - fnv (11) 
"^tiu = U^-h^u (12) 

where M"^ = + Mj + Ml and obviously + + ^m'' = 0' 

one can transform the kinetic term in (9) to 
which is not diagonal due to the last two terms. Before making it diagonal let us use 



^1 MgMfMh^ 



MMfh MMgf 
To make them diagonal one may now use the following transformations 



1 



to get 



M 



ib" = - {<b' \ = - lib' 



which is not yet canonical due to different weights. Eventually by assuming 



the kinetic term will be canonical as follows 

/iJ/'^ -■- ap T J-J-jui/'^ ■'"la/3 "nv(^ "apt 
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with new fields T, 11 and ^. Let us now use the above procedure for the potential term in Lagrangian (9). The 
corresponding potential term with a straightforward algebraic calculation is 

L (t\ Y\i^'' — TT'' TT'^ "1 - (" "^^^ _ \ — — (y\ "I^^ _ TT'' "1 



where 



Ml 



Ml 



M2 



M%M^f{ml + mi) + {Mf^ + Mgff{mlM^„) 
M%Ml^{ml + mi) + (M/„ - Mgf f{mlMl^) 



M'^M'jMl 



The full Lagrangian for these new variables is therefore 



(14) 



yyf 2 ti/t2 \ 



It is obvious that one of the gravitons, T, is masslcss though the other gravitons not only are massive but also 
interactive with arbitrary coefficients. This interaction term is a new prediction for having more than two gravitons 
and cannot be resolved by a new redefinition of the fields while taking the kinetic term canonical. This fact can be 
emphasized by noting that in the primary Lagrangian (9), three dimensional parameters exist. Consequently, these 
three dimensional parameters appear naturally after the field redefinitions. However, now two of them are represented 
as mass terms and the other as an interaction term. However, it is possible to resolve the interaction term by fine- 



tuning as m\ 



m\ and M\ 



g 



which means that in the primary Lagrangian (9) there are just two dimensional 



parameters. An interesting point is that the massless graviton corresponds to 



which is exactly the average of the primary fields in (9). To have an idea on the masses, let us assume the special 



case of rn? 



mt 



m| = m\ and ^/2>M = Mg = Mf = which means we have just one dimensional scale with 



the primary Lagrangian (4) being totally symmetric under changing the metrics. In this case the final gravitons 

have masses and = M| = 3 x ^ where ^ appears in (9) as mass tcrm[21]. The interesting point is that 
the coefficient "3" is exactly the number of metrics which exist in the model. This argument can be generalized for 
A/'-metric formalism by a little linear algebra for matrices. So in A'-metric formalism, one graviton is massless and 
other — 1 gravitons have masses as Mf = N x m?,yi e {1,2,..., AT — 1} where m? is the mass before the field 
redefinitions. 

It should be notic;cd that some of the above conclusion is true for the A^-metric formalism. In fact in A-metric 
models with ghost-free potential terms such as those in (9), after the redefinition of fields, the model represents a 
massless graviton which is the average of the metrics and the N — 1 massive gravitons with interaction. For a quick 
check let us innumerate the number of dimensional parameters. For the primary Lagrangian the number of interaction 
terms are (2)1 but after field redefinition one field becomes masslcss. The other A^ — 1 remaining fields are massive 
and present A^ — 1 dimensional parameters as well as ^) new interactive terms. So finally wc have A — 1 + ^) 
dimensional parameters which is exactly (^). In the language of [17], this model can be understood as follows; the 

primary Lagrangian (4) represents N sites which are linked to each other with links. This representation can be 
reduced to a representation with A sites with one site which has no link to the remaining A' — 1 sites. But these 
N — 1 sites are linked not only to each other but also have self- links or respectively have interactions and masses. 



A. Coupling to matter 



As mentioned in [7] , the coupling of multi- metric theories in the context of massive gravity to matter fields has not 
yet been solved. However, by considering the symmetry it seems as if one can constrain the coupling of the matter 
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to gravitons. By looking at (9) and assuming that there is no preference between the primary metrics one may then 
conclude that each kind of coupling is natural to be symmetric under the interchange of the metrics. In addition, from 
the linearized theory, it is obvious that there is a natural choice for the symmetric form of the metric which can be 
constructed out of their average. So, just by considering the symmetry, it seems that any correct coupling of matter 
should be to a symmetric combination of metrics which reduces to their average in the linear form. This argument 
becomes more significant if we note that the average of the metrics in the linear theory is simply the massless graviton. 
This means that whenever the linear theory is applicable, the matter couples to massless gravitons directly and does 
not see the massive ones. The existence of massive gravitons may show themselves in higher order terms which in turn 
may be verified in the context of cosmological perturbations. In the power spectrum of cosmological perturbations 
(e.g. curvature perturbation) there is no difference between Einstein gravity and massive gravity. However, in the 
bispectrum of perturbations (i.e. non-Gaussianity) the massive gravity should have different predictions. 

IV. COSMOLOGICAL SOLUTIONS 
A. Cosmological equations 

In this section we obtain the cosmological solutions of the above trimctric theory. In the massive gravity case, it 
has been shown that the flat FRW space-time is not the solution of the model [13]. However, open FRW solutions do 
exist [14]. Recently, the cosmological solutions for the non-linear ghost-free bimetric action was obtained in [15, 16]. 
Now we generalize them to the trimetric case. 

Let us assume that all three metrics can be described in terms of an isotropic and homogeneous space-times 

dsl = -N{tfdf + a{tf (dr^ + r^dif) 
ds} = -M{tfdt^ + b{tf {dr'' + r'^dfl'^) 

dsl = -Q{tfdt^ + w{tf {dr"^ + r'^dO?) . (15) 
Plugging these into action (4), we can read the reduced Lagrangian as follows 

Crea = 6 (^M|— +M]—+ Ml-^^ (16) 

+ 2mlM^f (^PoNa^ + j3i{Ma^ + 3Na^b) + 3j32ab{Ma + Nb) + j33b^{3Ma + Nb) + P^Mb^ 
+ 2mlM^^ (^^oNa^ + ^i{Qa? + SNa'^w) + ^^2aw{Qa + Nw) + •jsw^iSQa + Nw) + ^^Qw^ 

+ 2m|M|^ (^oMb^ + ai{Qb^ + 3Mb'^w) + 3a2bw{Qb + Mw) + azw^{3Qb + Mw) + aiQ'h 

The equations of motion can be obtained by varying the reduced Lagrangian with respect to scale factors with the 
result 




m J 

2m\Mlf (s^qNo^ + 3/3ia(Ma -|- 2Nb) + 3/326(2Ma + Nb) + d.^^Mb'^ 
2mlM^^ (^6joNa^ + 3jia{Qa + 2Nw) + '6j2w{2Qa + Nw) + ijsQw^^ = 0, (17) 



n.rl bbM 

6M? TT +2— -2 



- 2mlM^f (iPiNa^ + 3(32a{Ma + 2Nb) + 3P3b{2Ma + Nb) + S/SiMb^ 

- 2mlM]h (saoMb^ + 3aib{Qb + 2Mw) + 3a2w{2Qb + Mw) + 3a3Qw'^^ = 0, (18) 
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WW wwQ 



g2 J 

- 2mlM^f^ ^371 TVa^ + 3^2a{Qa + 2Nw) + 373u;(2Qa + Nw) + S'JaQw^ 

- 2mlM'jf^ i 3aiM62 + 3a2b{Qb + 2Mw) + 3a3w{2Qb + Mw) + ia^Qw"^ 



0, 



and variation of the Lagrangian with respect to the lapse functions become 



- 2m^M^^ ( 7oa^ + 371 wa^ + 3720^^ + 73^^ ) = 0, 



0, 



„ , ,9 WW 



2m^Mff^ ao6 + 3aiwb^ + 3a2bw^ + asw' 



- 2mlM^f^ ( 7ia^ + 372^0^ + 373aM;^ + 7411;' 



2mlMjf^ (^aib^ + 3a2wb'^ + Sa^bw'^ + a4W^^ = 0. 



(19) 



(20) 



(21) 



(22) 



B. Solution I 



In this section we are going to show that there is a solution for the above equations. To do this we suppose 

a = b = w (23) 

and N = M = Q = 1, that is, we work in the comoving gauge. Equation (20) reduces to two disjoint [22] equations 

d = ±uja, (24) 

1 

1 ■ ■ 



mfM^j ( 00 + 3^1 + 3/32 + /33 ) + m^M^\ ( 70 + 371 + 372 + 73 



and a ~ e+"* for minus and plus signs respectively. The latter, which is more 



Obviously, the solutions arc a 

interesting, is a de-Sitter solution with a cosmological constant A = 3w^. Now it should be checked if this solution 
also satisfies the first equation in (17) for a = 6 in the comoving gauge. This equation is 

d^+2aa = 3a;^a^ (25) 

where lo has been defined previously. Obviously, both of a = e~"* and a = e~^"* satisfy the above equation. Now we 
should check what are the conditions imposed on the remaining equations. One can simply find that the following 
equation should be satisfied due to (21) and (22) 



1 



V3M„ 



VSMf 
1 



mjM^f i^Po + 3^1 + 3/32 + Paj + mlM^h [lo + 37i + 372 + 73 
mlM^j (j3i+ 3/32 + 3/33 + /34 ) + ^3^/^ ( + 3ai + 3a2 + 03 



VSMh 



^2Mgh ( 71 + 372 + 373 + 74 I + rn^Mff^ ( ai + 3a2 + 3a3 + a4 



where the first line is the definition of uj and the last lines should be seen as constraints on the parameters such that 
a = b = w = e^"*. It is obvious that equations (18) and (19) are automatically satisfied as equation (17). 
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C. Solution II 



Let us now see if in such a massive cosmology with more than two metrics it is possible to have a solution with one 
Minkowski metric. This solution is impossible for the bimctric casc[23]. This assumption may be interesting because 
of its resemblance to ordinary massive cosmology [13] where a static metric together with a dynamical one exist. 

Using equation (22) and assuming that w is a constant, one gets the following relation 

which should be satisfied for all a and h. Since a and b are dynamical fields, the above condition is satisfied if and 
only if RHS and LHS are equal term by term. The general solution for this is as follows [24] 



a = 6, 



w = 1, 



mlM^h^i = -mlM%ai z e {1, 2, 3, 4}. 



(26) 



What is done here represents a possible cosmological solution for the trimetric formalism of massive gravity. This 

solution is deduced by assuming that one of the metrics is static. This property makes this solution interesting and 
comparable to the assumptions in [13]. Now, equation (20) with the above assumption becomes 



SM^ad^ = 2m\Mlf (^po + 3/3i + 3/32 + Pa^ + '^mlM^^ ^700^ + 3710^ + 372a + 73) , 
and equation (21) with a = b 

GMjad'^ = 2mlM^f ^/3i + 3^2 + 3^3 + + 2mlM'^^ (^aoa^ + Saio^ + 3a2a + as^ . 



(27) 



(28) 



By comparing the above two equations and our assumption of taking one of the metrics as static (26), it is easy to 
see that the following conditions have to be satisfied by the parameters 



1 

Ml 



mfM^f {^0 + 3/3i + 3/32 + /Sa) + mlM^hlo 

'yi = ai = 0, 



Mj 



mlMlf (/3i + 3/32 + 3/33 + Pi) + mlM%ao 



ie {1,2,3}. 



Actually, with the above conditions one should consider the remaining equations to get a correct solution. However, 
there is no need to pursue the matter any further in this case since, as we saw, the necessary condition of taking one 
of the metrics static is the same as assuming 7^ = = for i e {1, 2, 3}. By looking at Lagrangian (16), it is easy 
to find that this condition prevents the third metric in (15) interacting with others. So this solution is the trivial one 
and is without any importance. 

Finally, we showed that as a result of [13], there is no non-trivial solution for the massive cosmology studied here 
if, at least, one of the metrics is static. 



V. CONCLUSIONS AND QUESTIONS 

We have considered the generalization of the work done in [7], from a bimetric theory to a trimetric theory. This 

trimetric theory in fact has many of the consequences of the TV-metric theory, with the advantage of simplicity. What 
we get in this paper, is the conclusion that, in the multi-metric theories, only one graviton becomes massless at the 
linear level, which is the average of the N metrics. The other N — 1 spin-2 fields become massive and in general 
interact with each other. This has a general conclusion that ordinary matter must couple to metrics in such a way 
that at the linear level only the coupling to the average metric survives. This constraint on the form of the metric 
may be important when one addresses the subtleties of the multi-metric formalism. In this formalism, the gcionictrical 
interpretation of metrics e.g. the meaning of the distance and covariant derivative (parallel transportation) are 
yet unresolved. But with the above restriction which comes from coupling to matter a large class of multi-metric 
formulations are ruled out. 

The cosmology of this multi- metric theory is also interesting and generalizes the result in massive cosmology [13]. 
We have seen in this paper that, if we take one of the metrics as Minkowskian, it will become impossible to get a 
non-trivial cosmological solution for the theory. This result is however true only in a flat case. It is a matter of further 
investigation to see if the theory can have non-trivial non-flat cosmological solutions, as in the massive cosmology 
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case [14]. It is also crucial to consider the cosmological perturbation for FRW cosmology in this context due to the 
existence of a lot of observational data. From a theoretical viewpoint, when having more than one metric, it is natural 
to get not only the adiabatic mode but also entropy perturbations without any need for extra matter fields. 

We should also mention that the multi-metric theory constructed here is not proven to be ghost- free [25]. In fact, it 
is possible that the interaction terms introduced here may not suffice to cancel all Boulware-Deser ghosts of all spin-2 
fields. This would require further investigation in a future work. 

Acknowledgments 

We would like to thank B. Bassett, P. Creminelli, C. de Rham and P. Moyassari for their comments and discussions. 
N. K. would like to thank P. Creminelli for very fruitful discussions during his stay at ICTP and the organizers of the 
Workshop on Infrared Modifications of Gravity. 



[1] M. Fierz and W. Pauli, On Relativistic Wave Equations for Particles of Arbitrary Spin in an Electromagnetic Field, Proc. 

R. Soc. London A 173, 211 (1939). 
[2] H. van Dam and M. J. G. Veltman, Massive and mass-less Yang-Mills and gravitational fields, Nucl. Phys. B 22, 397 

(1970); V. I. Zakharov, Linearized graviton theory and the graviton mass, JETP Lett. 12, 312 (1970). 
[3] A. I. Vainshtcin, To the problem of nonvanishing gravitation mass, Phys. Lett. B 39, 393 (1972). 
[4] D. G. Boulware and S. Deser, Can Gravitation Have a Finite Range?, Phys. Rev. D 6, 3368 (1972). 

[5] C. de Rham and G. Gabadadze, Generalization of the Fierz-Pauli Action, Phys. Rev. D 82 (2010) 4, arXiv:1007.0443v2 
[hep-thj. 

[6] K. Hinterbichler, Theoretical Aspects of Massive Gravity, arXiv: 1105.3735. 

[7] S. F. Hassan and R. A. Rosen, Bimetric Gravity from Ghost-free Massive Gravity, arXiv:1109.3515vl [hop-th], 

S. F. Hassan and R. A. Rosen, Confirmation of the Secondary Constraint and Absence of Ghost m Massive Gravity and 
Bimetric Gravity, arXiv:1111.2070vl [hep-th]. 
[8] N. Rosen, General Relativity and Flat Space. I, Phys. Rev. 57 (1940). 
[9] N. Rosen, General Relativity and Flat Space. II, Phys. Rev. 57 (1940). 
[10] I. T. Drummond, Bi-metric Gravity and "Dark Matter", Phys. Rev. D 63, 043503 (2001), arXiv: 0008234v2 [astro-ph]. 
[11] M. Milgrom, Bimetric MONO gravity, Phys. Rev. D 80, 123536 (2009), arXiv: 0912. 0790v2 [gr-qc]. 

[12] C. de Rham, G. Gabadadze and A. J. Tolley, Resummation of Massive Gravity, Phys. Rev. Lett. 106 (2011) 231101, 
arXiv;1011.1232 [hcp-th]. 

[13] G. DAmico, C. do Rham, S. Dubovsky, G. Gabadadze, D. Pirtskhalava and A.J. Tolley, Massive Cosmologies, 
arXiv:1108.5231vl [hep-th], 

K. Koyama, G. Niz and G. Tasinato, Analytic solutions in non-linear massive gravity, Phys. Rev. Lett. 107 (2011) 131101, 

arXiv:1103.4708 [hop-th], 

K. Koyama, G. Niz and G. Tasinato, Strong interactions and exact solutions m non-linear massive gravity, Phys. Rev. D 
84 (2011) 064033, arXiv:1104.2143 [hep-th], 

K. Koyama, G. Niz and G. Tasinato, The Self-Accelerating Universe with Vectors in Massive Gravity, arXiv: 11 10.2618 
[hep-th]. 

[14] A. Gumrukcuoglu, C. Lin and S. Mukohyama, Open FRW universes and self- acceleration from nonlinear massive gravity, 

JCAP 11 (2011) 030, arXiv:1109.3845vl [hep-tii]. 
[15] M. S. Volkov, Cosmological solutions with massive gravitons in a bigravity theory, arXiv: 1110.6153vl [hop-th], 

D. Comelli, M. Crisostomi, F. Nesti and L. Pilo, FRW Cosmology in Ghost Free Massive Gravity, arXiv:1111.1983vl 

[hep-th]. 

[16] M. von Strauss, A. Schmidt-May, J. Enander, E. Mortsell and S. F. Hassan, Cosmological Solutions in Bimetric Gravity 

and Their Observational Tests, axXiv:1111.1655vl [gr-qc]. 
[17] N. Arkani-Hamed and M. D. Schwartz, Discrete Gravitational Dimensions, Phys. Rev. D 69, 104001 (2004), arXivihep- 

th/0302110. 

[18] N. Arkani-Hamed, H. Georgi and M. D. Schwartz, Effective Field Theory for Massive Gravitons and Gravity in Theory 

Space, Ann. Phys. 305 96-118 (2003), arXiv:hep-th/0210184. 
[19] S. F. Hassan and R. A. Rosen, On Non-Linear Actions for Massive Gravity, arXiv: 1103.6055v3 [hep-th]. 
[20] Note that our definition of perturbations is difi^erent with [7] due to a factor. It does not affect the final results. 
[21] Note that the dividing by "2" is a consequence of our assumptions for this special case i.e. = m\ = m| = m| and 

VSM ^ Mg = Mf = Mh and definitions (5) 
[22] The hncar combination of their solutions is not a solution. 

[23] This is obvious from recent works on this topic [15] and [16] though they did not mention it. The procedure is similar to 
what we have done for trimetric model. 



10 



[24] One can assume a = C^b and w — c for constant and c. But these assumptions will reduce to (26) by a re-definition of a, 

and 7i Vi € {1, 2, 3, 4} without any physical consequences. 
[25] However it seems the mechanism for bimetric model (see the second reference in [7]) may be applicable. 



